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Abstract 

Production of J/i/j mesons in heavy ion collisions is considered within the 
statistical coalescence model. The model is in agreement with the experi- 
mental data of the NA50 Collaboration for Pb+Pb collisions at 158 A-GeV 
in a wide centrality range, including the so called "anomalous" suppression 
domain. The model description of the J/ip data requires, however, strong 
enhancement of the open charm production in central Pb+Pb collisions. This 
model prediction may be checked in the future SPS runs. 
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Typeset using REVTgX 



Production of charmonium states J / if) and if)' in nucleus- nucleus collisions has been stud- 
ied at CERN SPS over the previous 15 years by the NA38 and NA50 Collaborations. This 
experimental program was mainly motivated by the suggestion [|TJ to use the J /if) as a probe 
of the state of matter created at the early stage of the collision. The original picture |J (see 
also 0] for a modern review) assumes that charmonia are created exclusively at the initial 
stage of the reaction in primary nucleon-nucleon collisions. During the subsequent evolution 
of the system, the number of hidden charm mesons is reduced because of: (a) absorption 
of pre-resonance charmonium states by nuclear nucleons (normal nuclear suppression), (b) 
interactions of charmonia with secondary hadrons (comovers), (c) dissociation of cc bound 
states in deconfmed medium (anomalous suppression). It was found [[| that J J if) suppres- 
sion with respect to Drell-Yan muon pairs measured in proton-nucleus and nucleus-nucleus 
collisions with light projectiles can be explained by the so called "normal" (due to sweeping 
nucleons) nuclear suppression alone. In contrast, the NA50 experiment with a heavy projec- 
tile and target (Pb+Pb) revealed essentially stronger J /if) suppression for central collisions 
0-0. This anomalous J /if) suppression was attributed to formation of quark-gluon plasma 
(QGP) 0|, but a comover scenario cannot be excluded ||. 

A completely different picture of charmonium production was developed recently within 
several model approaches PHn|. In contrast to the standard approach, hidden charm mesons 



are supposed to be created at the hadronization stage of the reaction due to coalescence of 
c and c quarks created earlier. In this case the J /if) yield is not restricted from above 
by the normal nuclear suppression curve. Therefore, neither anomalous suppression nor 
enhancement are excluded. 

In the present letter we consider the statistical coalescence model (SCM) fllQ| , [TT|j of char- 



monium production. We assume that c and c are created at the initial stage of the reaction 
in primary hard parton collisions. We neglect creation of cc pairs after the hard initial stage 
as well as their possible annihilation. Then, the number of charmed quark- ant iquark pairs 
remain approximately unchanged during subsequent stages. They are distributed over final 
hadron states at the hadronization stage in accord with laws of statistical mechanics. The 
SCM provides an excellent quantitative description of the NA50 data on centrality depen- 
dence of J /if} production in Pb+Pb collisions at SPS, provided that the number of nucleon 
participants is not too small (N p > 100). The peripheral collision data can be explained 
qualitatively. 

If creation of heavy quarks is indeed a hard process only, the average number (cc)ab(&) of 
produced cc pairs must be proportional to the number of primary nucleon-nucleon collisions. 
Then the centrality dependence of (cc)ab(&) can be calculated in Glauber's approach: 

(cc) AB{b) = or» N T AB (b). (1) 

Here b is the impact parameter, 7ab(&) is the nuclear overlap function (see Appendix) 
and a^ N is the cc production cross section for nucleon-nucleon collisions. As discussed in 



Ref. deconfmed medium can substantially modify charm production in hard collisions 
at SPS. Therefore, a^ 5 N in A+B collisions can be different from the corresponding cross 
section measured in a nucleon-nucleon collision experiment. The present analysis considers 
a^ 5 N as a free parameter. Its value is fixed by fitting the NA50 data. 

Event-by-event fluctuations of the number of cc pairs follow the binomial distribution, 
which can be safely approximated by the Poisson distribution because the probability to 
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produce a cc pair in a nucleon-nucleon collision is small: 



P k (b) = exp -(cc) 



>AB(b) 



[(CC) 



AB{b) 



k\ 



(2) 



where Pk(b) is the probability to produce k cc pairs in an A+B collision at impact parameter 
b. Assuming exact cc-number conservation during the evolution of the system, the SCM 
result for the average number of produced J/ ip per A+B collision is given by [|Tl[ 



(J/i>)AB(b) = (cc)aBQ>) 1 + (cc) A B(b) 



Here 
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is the total open charm thermal multiplicity. The sum runs over all known (anti) charmed 
particle species [[Uj. The total J/ip multiplicity includes the contribution of excited char- 
monium states decaying into J/ip: 



AI tot 



E R U)Nr 



(5) 



i=j/-0,xi>X2,V'' 



Here R(j) is the decay branching ratio of the charmonium j into J/i/j: R(J/tp) = 1, R{xi) ~ 
0.27, R(x2) ~ 0.1A and R(ip') ~ 0.54. The multiplicities iVj are found in the grand canonical 
ensemble formulation of the equilibrium hadron gas model: 



N 3 = Vn j( T ,VB) 



djV c^l T m 2 T ^ 
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(6) 



Here V and T are the volumef] and temperature of the HG system, rrij and dj denote, 
respectively, the masses and degeneracy factors of particles. K 2 is the modified Bessel 
function. The chemical potential fij of the particle species j in Eq. (§) is defined as 



bjHB + sj/is + Cjfi C 



(7) 



Here bj,Sj 



and Cj represent the baryon number, strangeness and charm of the particle j, 
respectively. The baryonic chemical potential \xb regulates the baryonic density. The 
strange \l$ and charm fic chemical potentials are found by requiring zero value for the 
total strangeness and charm in the system . In our consideration we neglect small effects of 
a non-zero electrical chemical potential. 

We assume that the chemical freeze-out occurs close to (or even coincide with) the 
hadronization (where charmonia supposedly formed). Therefore for the thermodynamic 
parameters T and /ig we use the chemical freeze-out values found |L7] by fitting the HG 
model to the hadron yield data in Pb+Pb collisions at SPS: 



1 We use ideal HG formulas and neglect excluded volume corrections. 



3 



T = 168 MeV, \i B = 266 MeV 



(8) 



Uncertainties in the freeze-out parameters exist due to time evolution of the system through 



the phase transition [O] and because of possible change of effective hadron masses in hot and 



dense hadron medium [ 19| . To check robustness of the predictions, an independent parameter 
set PU| is also used. It has been obtained by assuming strangeness and antistrangeness 



suppression by factor j s : 

T = 158 MeV, fx B = 238 MeV 7, = 0.79 . (9) 

The system is assumed to freeze-out chemically at some common volume. This is fixed 
by the condition of baryon number conservation: 

N p (b) = Vn B (T,fi B )=V J2 b jnj (T,fi B ), (10) 

j=N,N,A,A,... 

Here N p (b) is the number of participating nucleons at impact parameter b, n B (T,fi B ) is the 
net baryon density at the chemical freeze-out. The sum in Eq.(|l^) runs over all (anti)baryon 
species. The linear relation (|10"D between N p and V as well as constant values of T and \x B 
parameters are assumed for all collisions with different values of impact parameter b. 

Eq. (H) gives the total number of produced J ftps. They decay into ji + \T with the 
probability B$> = (5.88 ± 0.10)% Q. Only the fraction r] of fi + fi pairs that satisfies the 
kinematical conditions 

0<y<l, (11) 
-1/2 < cos6 < 1/2 (12) 

can be registered by the NA50 spectrometer. Here y stands for the rapidity of a /i 4 \i~ pair in 
the center-of-mass frame of colliding nuclei. 9 is the polar angle of the muon momentum in 
the rest frame of the pair. An estimate of 7] is impossible without detailed information about 
the hydrodynamic expansion of the system and the conditions at the thermal freeze-out. We 
shall therefore treat rj as one more free parameter. 

In the NA50 experiment the Drell-Yan muon pair multiplicity (either measured or calcu- 
lated from the minimum bias data) is used as a reference for the J/ip suppression pattern. 
Similarly to cc pairs, the number Drell-Yan pairs is proportional to the number of primary 
nucleon-nucleon collisions: 

{DY') AB{b) = a™T AB (b), (13) 

where the nucleon-nucleon production cross section of fi + fi Drell-Yan pairs. The 

prime means that the pairs should satisfy the kinematical conditions of the NA50 spectrom- 
eter (|i~l|) and (0). As the Drell-Yan cross section is isospin dependent, an average value is 
used: 

NN °DY' 

For the case of Pb+Pb collisions, A = B = 208 and o p ^ = 1.49 ± 0.13 fib §]. 
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Hence, the quantity to be studied is the ratio 



H W - — Tnv^\ ~ V b w~nn \} + a cc i ab{o) , . , 2 • (15) 

\L>Y )AB(b) a DY' l iV 0/^J 

It is convenient to rewrite the last expression in a simpler form 

m = c l±^M (16) 

and treat C and a^ 5 N as free parameters. In this form our fitting procedure does not depend 

on chemical freeze-out conditions. The new free parameter C is connected to 77 by the 
expression 



Here we have introduced the total open (anti)charm density: n Q = N Q /V and total J /if) 
"density": rijL = N^h/V. The relation between C and r\ does depend on freeze-out 
conditions, but our calculations with the parameter sets flU) and @ have shown that this 
dependence is not essential. 

In the NA50 experiment, the neutral transverse energy of produced particles Et was used 
to measure centrality of the collisions. This variable, however, provides a reliable measure 
of the centrality only if it does not exceed a certain maximum value: Et ^100 GeV (see 



also Ref. pljp2| )■ To show this we have calculated the dependence of the average number 
of participants on the transverse energy N p (Et). 

The conditional probability to measure some value of Et at fixed impact parameter b is 
given by a gaussian distribution: 

P {ET \ b) , 1 ~ T f gL^wn (18) 

j2*q'aN r (b) \ 2q 2 aN r (b) j 

Analyzing the experimental situation, we are interested in a quite opposite question: how 
events with fixed Et are distributed with respect to the centrality. The answer is 

pmF , _ bP(E T \b)P„„(b) m] 

where P int (b) stands for the probability (see Appendix) that two nuclei at fixed impact 
parameter b interact (at least one pair of nucleons collides). The average number of partic- 
ipating nucleons at fixed Et is then given by the expression: 

^■■r^'- ^M 1 (2o) 



The parameter values q = 0.274 GeV and a = 1.27 |23[ are fixed from the minimum bias 
transverse energy distribution. 
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The result is shown in FigfJ. As is seen, the transverse energy is simply related to the 
number of participants Et = qN p in the domain Et ^ 100 GeV. Outside of this domain 
N p does not change essentially as E T grows. Therefore the data at E T > 100 GeV do not 
represent centrality dependence of the J/i[) suppression pattern but rather its dependence 
on fluctuations of the stopping energy at fixed number of participants. In principle, influence 
of such fluctuations on J j ip multiplicity can be studied in the framework of our model, but 
information concerning the corresponding fluctuations of the chemical freeze-out parameters 
T and hb would be needed. Experimental data that would allow to extract this information 
(hadron yields at extremely large transverse energy) are not available at present. Therefore, 
we restrict our analysis to centrality dependence of J / ip production and do not use the data 
corresponding to large transverse energies Et > 100 GeV. 

On the other hand, the SCM is not expected to describe small systems. This can be seen 
from ip' data flO| . In the framework of SCM the multiplicity of ip' is given by the formula @ 
with the replacement Nj°L — > N^> . Therefore, the ip' to J j ip ratio as a function of centrality 
should be constant and equal to its thermal equilibrium value. The experimental data |24[ 



(see also a compilation in Ref. |10[) are consistent with this picture only at rather large 
(N p ^ 100) numbers of participants p5| . 

Hence, the applicability domain of the model is limited to 

27 < E T < 100 GeV. (21) 

Note that the most precise and abundant NA50 data (see Fig. [T]) correspond to this kine- 
matical region. 

At Et ^ 100 GeV the formula ([16]) and the equation 

E T = qN p (b) (22) 

give a parametric dependence of the ratio R on the transverse energy. This dependence for 
the parameter set 

C = (2.59 ±0.25) • 10 3 
ag N = (34 ± 10) fib (23) 

is plotted in Fig. [j]. The free parameters were fixed by fitting three sets of NA50 data ||[7| 
within the applicability domain fl2~T|) of the model by the least square method. The model 
demonstrates excellent agreement with the fitted data (x 2 /dof = 1.2). 



Extrapolation of the fit to peripheral collisions reveals sharp increase of the ratio (15) 
with decreasing N p . Such behavior in the SCM can be understood as the following. The 
smaller is the volume of the system the larger is the probability that c and c meet each 
other at hadronization stage and form a hidden charm meson. As is seen from Fig. [I], 
this is not supported by the data: the SCM curve lies above the experimental points in 
the low Et region. On the other hand, the normal nuclear suppression model also fails to 
explain the leftmost point from the 1996 standard analysis set and two leftmost points from 
the 1996 minimum bias set. Those theoretical calculations underestimate the experimental 
values. It is natural to assume that an intermediate situation takes placefl Some fraction of 



2 Similar combination of standard and SCM production mechanisms has been considered in Ref. 
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peripheral Pb+Pb collisions result in formation of deconflned medium. In these collisions 
charmonia are formed at the hadronization stage, and their multiplicities are given by SCM. 
The rest collisions (we shall call them 'normal collisions') do not lead to color deconfinement, 
therefore charmonia are formed exclusively at the initial stage and then suffer normal nuclear 
suppression. The experiment measures the average value, which lies between the two curves. 

The fraction of 'normal' events decreases with growing centrality. Their influence on J /if) 
production becomes negligible at N p ^ 100. To check this we repeated the above fitting 
procedure using only the experimental data corresponding to N p > 200. The quality of the 
fit is only slightly better: x 2 /dof = 1.1, the parameter values C = (2.73 ± 0.40) • 10 3 and 
a c5 N = (31 =t 12) fjh are consistent with the analysis of the full data set (|23|). 

Our picture is also supported by if)' data. The normal nuclear suppression influence 
nascent charmonia before the formation of meson states. Therefore its effect on if)' is the 
same as on J /if). The multiplicity ratio of if)' to J/if> in 'normal' nuclear-nuclear collisions 
should be the same as in nucleon-nucleon collisions and should not depend on the centrality. 
In the framework of SCM, the if)' to J/if> ratio, as was explained above, should be equal to 
its thermal equilibrium value, which is a few times smaller than the corresponding value for 
'normal' collisions. As the fraction of 'normal' events decreases, the measured ratio should 
decrease and then become constant and equal to its thermal value. The experimental data 
24[ indeed demonstrate such behavior [f^ , |10fl . 

The present analysis predicts strong enhancement of the total number of charm. From 
a pQCD fit of available data on charm production in p+N and p+A collisions, one could 
expect dcc N ~ 5.5 fib at y/s = 17.3 GeV. Our result (|23f) is larger by a factor of 4.5 -j- 8.0, 



which is around the upper bound of the charm enhancement estimated in Ref. |jT3 . 

Formation of deconfinement medium can change not only the total number of charmonia 
and open charm particles but also their rapidity distributions. For direct charmonium 
production in hard parton collisions, dimuon pairs satisfying the kinematical conditions ([TT|) 
and ( |T2"D account for a fraction of about r] hard « 0.24 in the total number of pairs originating 
from J /if> decays. (The value was found using Schuler's parameterization f27j] .) Our result 
(^3f) corresponds to r\ m 0.14, which is by a factor of about 0.6 smaller. This difference can 
be attributed to broadening of the J/ if) rapidity distribution. It is natural to expect similar 
modification of the open charm rapidity distribution. Because of this modification the open 
charm enhancement within a limited rapidity window can, in general, differ from the one for 
the total phase space. Assuming that the broadening for the open charm is approximately 
the same as that for J/if>, one obtains open charm enhancement by a factor of about 2.5 -=-4.5 



within the rapidity window fliTD , which is consistent with the indirect experimental result 



[26] for central Pb+Pb collisions. It was not, however, checked whether this approach is able to 



describe the centrality dependence of the J /if) suppression pattern in (semi) central collision region. 



3 Our previous study jn] were based on the J /if) multiplicity data f29| , which were extracted 
from the NA50 data |||?J assuming narrow rapidity distribution of J /if). In this case, the charm 
enhancement in the total phase space is by a factor of about 2-3.5 and does not differ from the 
enhancement in the limited rapidity domain, but the number of cc pairs should grow faster than 
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In conclusion, we have shown that the NA50 data on centrality dependence of the J /if) 
and if)' production in Pb+Pb collisions @,0,^4|] are consistent with the following scenario: 

The deconfined medium, which is formed in a Pb+Pb collision, prevents formation of 
charmonia at the initial stage of the reaction. Instead, hidden charm mesons are created 
at the hadronization stage due to coalescence of created earlier c and c quarks. Within 
this scenario, the color deconfmement does not necessary lead to suppression of J /if). Both 
suppression and enhancement are possible P0| . If the number of nucleon participants is not 
too small (N p > 100), the number of produced J ftp is smaller than in the case of normal 
nuclear suppression, therefore anomalous suppression is observed. As color deconfmement 
is present in most collision events for N p > 100, our model reveals excellent agreement with 
the experimental data in this centrality domain. The statistical coalescence model does not 
describe the NA50 data for the peripheral Pb+Pb collisions. It seems that the fraction of 
events producing the deconfmement medium is not dominating there and most of peripheral 
collisions follow the normal nuclear suppression scenario. Still, the presence a fraction of 
abnormal events could reveal itself in the deviation of the J/if> data up from the normal 
nuclear suppression curve. 

Our model analysis predicts rather strong enhancement of the open charm. This effect 
can also be related to the color deconfmement [TE\. The enhancement within the rapidity 
window < y < 1 is consistent with the indirect NA50 data [28|. A direct measurement of 
the open charm would be very important for checking the above scenario. 
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APPENDIX: NUCLEAR GEOMETRY 



The spherically symmetrical distribution of nucleons in the Pb-208 nucleus can be pa- 
rameterized by a two-parameter Fermi function (this parameterization is also known as 
the Woods-Saxon distribution): 

(Al) 

with c ~ 6.624 fm, a ~ 0.549 fm and po is fixed by the normalization condition: 



P(r) = Po 



1 + exp 



r — c 



the number of nucleon- nucleon collisions. Either data on rapidity distribution of J /if) or precise 
data on centrality dependence of the open charm multiplicity would help us to decide, which of 
these two versions is preferable. 
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47T 



drr 2 p{r) = 1. 



(A2) 



The nuclear thickness distribution T A (b) is given by the formula 

T A (b) = r dzp (VtfT~z 2 



and the nuclear overlap function is defined as 

Tab(P) = £^ dx J™ T A (y/^ + Z/ 2 ) T B ( + (i - b? 



(A3) 



(A4) 



From Eq.(|A2|), one can deduce that the above functions satisfy the following normalization 
conditions: 



/■oo roo 

2tt / dbbT A {b) = l, 2vr / dbbT AB {b) = 1. 
j o j o 



(A5) 



In Glauber's approach the average number of participants ('wounded nucleons') in A+B 
collisions at impact parameter b is given by 



NJb) = A 



dx 



+ B 



dx 



dy T A ^x 2 + y 2 j |l - 1 - a^T n ( V <V-' - (,, - />)'- 

~ 00 ^T B ( v /W(2/- & ) 2 m 



(A6) 



1 - a^T A ( Jx* + 



Here a™^ is the nucleon-nucleon total inelastic cross section. 

At large impact parameter, the nuclei may do not interact at all. Therefore N p (b) — > at 
b — ► oo. If one interested in the average number of participants, provided that an interaction 
between two nuclei has taken place, the relevant quantity is 



where 



N p (b)=N p (b)/P int (b), 



P int (b) = 1 - 1 - <nelT A B{b) 



AB 



(A7) 



(A8) 



is the probability for nuclei A and B to interact at impact parameter b. Although N p (b) 
differ from N p (b) at large b: N p (b) — > 2 at b — > oo, they are almost identical for more central 
collisions. 

The average number of nucleon-nucleon collisions can be calculated from 



N coU (b) = ABa^ el T AB (b) 



(A9) 



Provided that an interaction between two nuclei has taken place, the above formula should 
be modified as 



N coll {b) = N co u{b)/P int {b). 



(A10) 
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FIGURES 




FIG. 1. The dependence of J ftp over Drell-Yan multiplicity ratio on the transverse energy. The 
normal nuclear suppression curve is obtained at a a b s = 6.4 mb, where a a bs is the absorption cross 
section of preresonant charmonia by nuclear nucleons. Two vertical lines show the applicability 
domain of the model under consideration, see the text for details. 
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FIG. 2. The dependence of the average number of nucleon participants on the transverse energy. 
The dotted straight line corresponds to Eq.(f22]), 
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